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Abstract 

The Hubbard model is used to study an electronic system. In this paper we present 
the new path integral representation for Hubbard model. We have constructed the 
new supercoherent state which appears from a set of eigenfunctions of atomic limit of 
strongly correlated systems. Exact calculation of nonlinear representation of a super- 
group has been carried. This group defines the transformation of atomic base. The 
general formalism we elaborate for Hubbard model is the one widely used in the gauge 
field theory of the nonlinear representation of a superconformal group. 
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1 Introduction 



The Hubbard model was originally constructed to describe a metal- insulator transition 
for spin- dependent fermions in a simple way [U El E] . 

Today this model is still remain the main workspace for investigation of the strong elec- 
tronic correlation. There exist many approaches to this model for describing many electron 
system: the band limit approximation for the weak interaction between electrons and the 
atomic limit for electrons with the strong coulomb repulsion. We start the series of papers 
in which we are intending to elaborate the new approach to the Hubbard model. 

We will present a new path integral approarch to the strong interaction regime. The 
main ingredients for us is the usage of a supercoherent state acting upon an atomic base 
and the work with effective functional for an electronic system. We will develop the pro- 
cedure of geometric quantization for strongly correlated systems of electrons. Let us make 
the brief sketch of the program firstly developed in the series of works of one of the authors 
[I]. One of the main distinction of this approach to the Hubbard model is that we treat the 
local supergroup of the local space-time as the main object of our theory. This supergroup 
generates the transformation of the space-time coordinates which assign the arguments for 
any function describing this system. This representation of a supergroup in the superspinor 
space has to contain a Lorents or an SO (4) subgroup equal to the even subalgebra of the 
Hubbard operators. Odd Hubbard operators produce some superextension of the Lorents 
group into some supergroup. This dynamic supergroup is given by local dynamic superfields 
describing the local degrees of freedom of the strongly correlated electron system. We shall 
introduce the supercoherent state depending on generalized angles equal to the bosonic and 
fermionic fields of the system. Parametrised by x,y,z,t space-time manifold which deter- 
mines the arguments of wave function defines us the spinor and superspinor bundle. This 
superbundle is determined by the supercoherent state. 

There are two kinds of the gauge fields in superspinor bundle: one sort of fields is the com- 
posite fields equal to the quadratic combination of odd grassmann fields and other sort are 
nonlinear fields determining the local coordinates frame of the four-dimensional space-time. 
In general, the local superspinor bundle defines the nonlinear representation of a super- 
conformal group as a maximal group of 4-dimentional space-time of interacting fermionic 
system. 

As the first step of quantization of electronic system with strong coulomb repulsion 
we perform the reformulation of Hubbard model into the atomic limit formalism. This 
approach is well known in Mott-Hubbard insulators theory. We want to point out that our 
approach includes all elements of geometric quantisation [5]: for example possess the algebra 
of 4-dimensional rotations (Lorents) group, Cartan differencial one-forms which give us the 
lagrangian of system, the nonlinear representation of underlying supergroup as a ground for 
the supercoherent state. 

In this first paper of this series we will fulfill the exact calculation of the nonlinear 
representation of those supergroup generators of which appear in operator approach of atomic 
limit of the Hubbard model. This gives us the possibility to introduce the supercoherent 
state for the strongly correlated electron system and as a result to obtain the effective action 
in a future paper. 

Let us give the brief description of the calculation method for the finding representation 
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of dynamic supergroup in strongly interacting models. Our task is to find in this model such 
group structure which could help us to describe the specificity of strong correlation. We take 
the following construction as a base: 

1) we will collect space coordinates together with the time coordinate and will consider 
some curved space-time as a base in which Lorents subalgebra of superconformal group act 
in the spinor base on dynamical fields. 

2) full bases of electronic operators gives us the supergroup which is parametrised by 8 
dynamical fermionic fields: this fields comprise the conformal spinor. 

3) the superspinor representation of a supergroup gives us the supercoherent state de- 
scribed by the nonlinear function over odd grassmanian fields. This function characterises 
the local properties of the strongly correlated system. 

2 Atomic description of Hubbard model 

We consider the Hubbard model: 

ija i,a <j,i 

here a^a a j -electron creation and annihilation operators. n a ^ - electron density 
operator W, U, -band width, one- site electron repulsion and chemical potential. 

At first we represent the Hubbard model in atomic bases which determine the atomic 
limit. This limit appear as a result of the following procedure. A zero approximation of the 
atomic limit is described by one-site repulsion term: 

This hamiltonian can be diagonalized by the following one-site atomic eigenf unction: 

|0 |+ >-= <xf|0 |- >-= aj~|0 |2 >-= a+a+|0 y . (2) 

This bases gives us the fundamental representation of some supergroup in the space of 
dimension (2,2). Point out that the states |+ >-;|— are fermions but |0 >-; |2 >- 
are bosons in our construction. We insert later odd grassmann fields and make this fact 
obvious ( ie states |+ ]— will be depend on odd order of the grassmann fields but 
states |0 >~; |2 y -on even order of the grassmann fields). All operators in this bases will 
be the matrix which are determined by commutation and anticommutation relation giving 
some superalgebra. Full set of the Hubbard operators have 16 operators part of which 

(X 0+ , x°~, x +0 , x~°, x +2 , x~\ x 2+ , x 2 -) 

are the fermionic operators, but other part 

(X + -,X~ + , X ++ - X", X 02 , X 20 , X 00 - X 22 ) 

-are the bosonic operators. X*- 7 -Hubbard operators contain only one non-zero element equal 
1 sitting on site in the matrix representation. Point out that this set of operators gives 
some bases for some superalgebra. 
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We have the following representation for creation-annihilating operators in this bases: 

a + = X +0 + X 2 -a+ = X~° + X 2+ (3) 
The Hubbard model in this representation has the form: 

H=Uj2x? P --Wj2 X rX? (4) 

i,p ija/3 

3 Supercoherent state for Hubbard model 

In constracting the supercoherent state we use the following interesting observation in in- 
terpretation of the set of atomic operators and function for on-site Hubbard repulsion. This 
observation can be formulated as the following statement: six even Hubbard operators con- 
stitute the subalgebra isomorphic with algera of Lorentz group or algebra of four dimensional 
rotation group in spinor representation. Complete derivation of this statement will be obtain 
in subsequent paper. 

To characterise the state of the system by coherent state we input some fields which 
depend on coordinates x, y, z and the time t. We have three component dymanic vector of 
the electrical field 

E = (E + (x, y, z, t), E~(x, y, zt), E z (x, y, z, t)), 
three component dynamical vector of the magnetic field 

h = (h + (x,y,z,t),h~(x,y,z,t),h z (x,y,z,t)) 

and four component dynamical odd grassmann fields 

Xl(x, y, z, t), x* 2 (x, y, z, t), x* 3 (x, y, z, t), xl(x, y, z, t), 

which are the fermionic fields givving the components of maiorana field. All dynamical fields 
appear in supercoherent state in the followng manner: 

| G >= exp 

| G > Exponent here act in space of atomic eigenfunction (| >,| + >,| — >,| 2 >), 
Function | > is higthest weitght vector of suprgroup which representation is given by 
exponetial. 

4 Evolution operator for electronic systems 

The transition amplitude of the evolution operator of the quantum systems is given by 
the following expresion: < Z/le -1 ^*^ - **^^ >. We want to obtain the expression for the 
effective functional using the states \Z >. Time evolution of the system is given by the 
following operator: 



E z 

X*2 



E + 
h z h+ 
hr -K 



0>, 



(5) 



E~ -Xl Xl ~E Z 
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U(t,t Q ) =T ord exp(-i H(r)dr) 

J t 

if t — to — St is small, ie St « 1, then 



rto+St 

U(t + St,t ) = l-i H(r)dT 

J t 



It is follow from this expression that the symbol for evolutionary operator has the fol- 
lowing form: 

rto+St 

U(Z,Z*\t + St,t ) = exp(-i / H(Z,Z*\T)dr). 

J t 

We devide time interval [to,t] by the number N and obtain N small intervals for finding 
the expression for symbol U(Z, Z*\t, t ). Consider the matrix elements of evolution operator 
exp(—iH(tf—ti)) between the states < Zf \ and \Z { >. Factorising operator exp(—iH(tf—ti)) 
by inserting the identity operator f d/i(Z)\Z >< Z\ — 1 we obtain the following representa- 
tion: 

N 

< Zf\exp(-iH(t f -U))\Zi >= / Y[d/i(Z k ) < Z f \Z N > 

J k=i 

< Z N \e-^ H \Z N ^ > .... < Z k _ x \e^ H \Z k > ... < Z 1 \e~ ieH \Z i >, 

here e = ^p-.In first order of e we can transform this formular and place the symbol of 
operator H in the exponent 

<Z k+l \e^\Z k > = <Z k+1 \(l-ieH)\Z k > = + ^ (g) 

< z k+ i\z k > < Z k+ i\Z k > 

As a result we obtain the representation 



rr, fJ-L -i f \- N <Z k + l\ H \ Z k> 

<Z f \e- lH ^-^\Zi>= lim Y[dfi(Z k )<Z k+1 \Z k >e ^ k=1 <W*> , (7) 

N—>oo I - LJ - 
J k=l 

here \Z >= \Zi >; < Z N+ i\ =< Zf\. Let define a variation of the following type \Z >: 
\SZ k+ i >= \Z k+ i > — \Z k > . We have: 

N 

< Zfle'^f-^lZi >= lim / TT[rf/x(Z fc ) < Z k \Z k >] < Z f \Z N > 

N—>oo I 

J k=l 



exp(y(Ln(l — — ) - ie — — ). 

^ < Z k \/j k > < z k \z k _i > 

in linear-slice approximation. We take the following expresion for the time derivative: 
Considering the first order in e we can take the folllowing expresion for the time derivative: 



4 



d\Z > _ \SZ > 
dt e 

In first order in e, we obtain the final path integral representation of the evolutionary 
operator in coherent state formalizm 



p\Z(tf)>=\Z f > 

<Z f \e- iH{t f- t ^\Z l >= / D(Z,Z*)e- tS ^ z ^; 

J\Z(ti)>=\Zi> 



S[Z, z* 



[ * dt [ dr(- 
Jti Jv 



< Z(r,t)\i£f-H\Z(r,t> 



(8) 



< Z(r,t\Z(r,t> 
-i[Ln(< Z f \Z(t f ) >) - Ln(< Z t \Z(t t ) >)]). 
Measure of integration is given by the following expresion: 

D[Z,Z*]= JJ l[diM[Z(r,ty,Z(r,t)]<Z(r,t)\Z(r,t)>. 

ti<t<t f r 

This form of path intergal representation will be the starting point of our consideration. 



5 Nonlinear representation of supergroup in Hubbard 
model 

In construction of supercoherent state we have the followng supermatrix which we must 
compute analyticaly: 



U = exp 



/ E z E+ \ 

Xi K h+ 

X2 hr -h z 

V E~ X 3 Xa -E z J 

In this expression we have the fields of different statistics: for example, set of the fields 



x , V, z ), X2(t, x, y, z), xs(t, x, y, z), x*(t, x, y, z)) 

are the odd grassman valued function of space-time coordinates and describe the fermionic 
degree of freedom but fields 



(E z (t, x, y, z), E + (t, x, y,z),E (t, x, y, z), h z (t, x, y, z), h + (t, x, y,z),h (t, x, y, z)) 

describe the electro-magnetic degree of freedom, equal to two three component vectors of 
the space-time coordinates and are bosonic. In this paper we concentrate in calculating exact 
representation of exponent of the supermatrix in the coherent state. We take the dynamic 
electrical, magnetic and grassmann fields which depend on coordinates of 4-dimensional 
space-time manifold on definite the coordinates and omit (t, x, y, z) coordinates in subsiquent 
formulas. 
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Our general strategy will be to expand the supermatrix to N-order in fields. Them we 
can isolate and collect certain series and get some recurrent formular for general term in 
infinite series. Using this formula we can sum all terms to anytical compact representation. 
Analytical representation of the supermatrix elements will be the final point of our work. 
As a starting point we have the following exponencial expression for the representation of 
super extension of the Lorentz group in the spinor representation. 

Expanding this exponent in series we can obtain first and second order in the parameters 
b and h. We see that the polynomial series on the grassmann numbers can be classified 
in grassmann order n. All the supermatrix elements can be represented as a coefficients in 
grassmann polynomials of order n, where n=0, 1,2,3,4. 



6 Matrix series for the nonlinear representation of 
supergroup 

First of all point out that the supermatrix in exponent have two submatrix: one is the odd 
grassmann matrix and the other is the even submatrix containing only the fields of type 
Ei(t, x, y, z) and hi(t, x, y, z) type. As a first step we make expansion of the exponent for 
the even matrix. Expanding in series this exponent we can obtain first, second order and 
3,4,5 order in bosonic fields Ei,i 
n=0, 1,2 has the following form: 



1,2,3. and hi,i = 1,2,3. For example, the series for 



\ 
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h A h + 

-h A h z 
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here we introduce the following abbriviation b = \/ E 2 + E + E~, h = \Jh 2 + h + h~ 

For the odd grassmann number we have following expansion series of exponent. We write 

here two terms for the grassmann fields (xi, Xii X3? Xi) f° r obtaining coefficients in higher 

order in E~ and h; . 
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/ \ 

Xi o 

X2 o 

V o -xs X4 o / 



+ 



1 

2! 



/ o 

Xi-E* + X2h + + 
X2# 2 + Xih~ - X2K 

V 



-X3E+ xiE+ \ 

XiE + 

x 2 ^ + 

XsK + X±h~ - Xzh z -XaE z - X3h + - X4^ z / 



For even order of the grassmann variables we have following series for the composite 
bosonic fields: 



/ \ 





\ x20 / 



+ 



3! 



/ £+x20 \ 

-XiXsE + X iX±E + 

-X2X$E + X2X±E + 

\ he E + x2Q J 



+ 



1 

4! 



/ E+he + E+E z x20 

E + heh -E+h + x20 

-E + h~x20 E+he 

\ x20b 2 + h 2 x20 + E-E+x20 



E +2 x20 \ 



E+he - E+E z x20 J 



here we introduce the following abbriviation: x20 = — XzXi + X4X2, = —h z X2Xi + 
fr~XiX4 - h z X2Xi ~ h + X2Xs, 

heh = -h z xiX3 + h'XiXi ~ h z XiX3 ~ h + X2X3 

The expansion series for third order of the grassmann variable has the following form of 
first two terms 



1 

4! 







V 



+ 



\ 

-al3x 4 £ + 

-al3x 3 £ + 

a33x 2 £ + a33xi£ + / 

-a33x 2 £ +2 -a33 X iE+ 2 

-al3xiE + E z -al3x 4 £ +2 

-al3xsE + E z -al3x 3 E +2 

-a33x 2 E + E z ~a33xiE + E z 

here al3 = X1X2, a33 = XzX*- 

Forth order of grassmann variable has the following series of three terms: 



1 
5! 



/ 



V 



1 

5! 



/ \ 

000' 



y 2c4hiE+ / 



1 

+ 1 



/ 2c4hiE +2 \ 

-cAhiE+ 2 

-c4hiE+ 2 

\ 2cAhiE+ 2 J 



+ 
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1 

7! 



f 2cAhiE +2 E z 2cAhiE +3 \ 

-cAhiE +2 h z -cAhiE +2 h+ 

-cAhiE +2 h~ cAME +2 h z 

\ AcAhiE + b 2 + 2cAhiE~E+ 2 + 2cAhiE+h 2 -2cAhiE +2 E z ) 



here cAhi = X1X2X3X4 

Expanding this series on 12 order in even and odd parameters we can obtain the analytical 
representation for the supermatrix elements. 

7 Analytical representation of supergroup 

Having series representation for supergroup for high order in fields we can obtain the general 
analytical form for the matrix elements. Our task is to obtain exact dependence of the 
matrix element over component of fields: E Zt E + ,E~;h z ,h + ,h~ but not dependence over b 
and h. 

For example the general form of Uu supermatrix elements in orders higher than 8 as a 
functions of the dynamic even and odd grassmanian fields are given by following expressions 

Ull = a° n + E z b° n + £ + (-XsXi + X4X 2 )K1 + E z b\\)+ 
E + (-h z X3Xi + ^~X4Xi - h + XsX2 - ^X4X2)(an + #*&n)+ 

XlX2X3X4(£ + ) 2 K 1 +£ 2 ^ 1 ) 

here the coefficients a\ and b\ are some series in b and h variables: 
a^, b®! have some infinite series over b and h variables. 

Let introduce following abbriviation for coefficients which appear in the suprmatrix ele- 
ments: 

a\\ = J2EE- b\\ = flEE; a 22 = flEE; b 22 = fEE; a\ x = DEDEfl; b* n = 
DEDEf 

Another matrix element equal to 

u 21 = (a° 21 + E z b° 21 ) X i + (h zX i + h + X 2){a\ l + E z b l 21 )- 

E + XiX2xMli + E z b 3 21 ) 
and we have some function for the coefficients: 

a° 21 = fA; a° 21 = fA; b° 21 = /3; a' 21 = /3; b\, = /2; a 3 21 = flEE; b 3 21 = fEE 
For u 3 i we have 

«3i = + E z b° 31 ) X 2 + (~h zX2 + h~xi){a\i + E z b l 31 )- 
E + XiX2xMi + E z b\ l ) 
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and for coefficients: 

a§i=/4; 6§i = A <4 = -/3; & = /2 
Last element of first column is 

«4i = + (-X3Xi + X4X2)(a 2 4 { + E + E-bl\)+ 

(-h z X3Xi + h~x±Xi ~ h + XsX2 - h z x±X2){af 1 + E + E'b 22 1 ) + 

2 X iX2X3X*E + (ai 1 + E+E'biJ 
For coefficients in this expression we have: 

b° u = sinh(6)/6; a\\ = /3; 6ft = flEE; 

a 22 = A &2? = al = fEE; b\ x = DEDEf 

For second column we have 

U12 = -E + a° 12 X3 + E + (h z X3 ~ h~x^)a{ 2 - 
(E + ) 2 X 2X3X4a 3 12 

For coefficients a\ and b\ : 

a? 2 = A a\ 2 = -/2; a? 2 = /EE 

For -u 2 2 we have: 

M22 = a° 2 + h z b° 22 - E + xi X 3a 2 22 + E + h + h~ (-XsXi + X^ajl 
E + (-h z XiX3 + h~XiX4 - h + X2X3 ~ KXiX^H + h z b 2 22 )- 
XiX2X3X4(E + ) 2 (a A 22 + h z b\ 2 ) 

For coefficients: 

a\ 2 = -f1- af 2 = fhh; af 2 = flhh; b% = fhh; a\ 2 = /3; b\ 2 = /2 
For -u 32 we have: 

m 32 = h~bl 2 - E + X2X3a 2 3 l + E + h~(-X3Xi + X±X2){af 2 + h z bf 2 ) 

E + h~(-h z X2X4 + h~xiX4 ~ h + X2X3 ~ KX2xMf2~ 
X1X2X3X A(E + ) 2 h~ a\ 2 

For coefficients: 

af 2 = f2; al% = flhh; bf 2 = fhh; af 2 = fhh; aj 2 = DEDhf2 
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For U42 we have: 

U42 = (a° 2 + E z b° 42 )x3 + (h z X3 - h~X4)(a\ 2 + E z b\ 2 ) + 

E + X2X3X4(at 2 + E z b\ 2 ) 

For coefficients: 

= b° 42 = f3; a\ 2 = /3; b\ 2 = /2; < = 6| 2 = /££ 

For -u 13 we have: 

«i 3 = E + a° n X4 ~ E + (h z X4 + h~X3)al 3 - 
(E + ) 2 Xi X 3X4a 3 13 

For coefficients: 

a? 3 = /3; a} 3 = /2; a? 3 = fEE; a\ 3 = fEE 
For -u 32 we have: 

U32 = h+b° 22 + E + XiX4a 2 2 \ + E + h + (- X 3Xi + X4X 2 )(a| + h z bf 3 ) 

E + h + (-h z xiXs + h~XiX4 - h + X2Xs ~ KxxXz)^- 
XiX2X3X4(E + ) 2 h + al 3 

For coefficients: 

a| = /2; a 2 i = flhh; b 2 £ = fhh; a\ 3 = fhh; a\ 3 = DEDhfl 
For M33 we have: 

U33 = a°22 ~ h zb°22 + E+ X2X4a 2 33 + E + h~h + (-X3Xi + X4X2V33 
E + (-h z X2X4 + h~xiX4 - h + X2X3 ~ KX2X4){a 2 3 2 3 - h z b 2 33 )- 
X\X2X3X4[E + f[a\ 3 - h z b\ 3 ) 

For coefficients: 

af 3 = f2; a\ 3 = fhh; a 22 = Dfhh; b 22 = fhh; 
a\ 3 = DEDhf3; b A 33 = DEDhf2 

For M43 we have: 

U43 = (a% + E z b° A3 )X4 + (h z X4 + h + X3)(a\ 3 + E z b\ 3 )- 

E + XiXzX4{a\ 3 + E z b%) 
10 



For coefficients: 

a% = f±, b° 43 = /3; a\ 3 = /3; ^3 = /2; <& = b\ 3 = -fEE 
For u 14 we have: 

u 14 = + (^ + ) 2 (-X3Xi + X4X 2 )a?4 + XiX2X3X4(£ + ) 3 < 

For coefficients: 

afi = flEE; a\ 4 = DEDEf 

For -u 2 4 we have: 

« 24 = E+b° 24 Xi + (h z Xi + h + X2)E + b l 2i - 
(E + ) 2 XlX 2X4all 

For coefficents: 

a° 24 = f2; b° 24 = f3; b\ 4 = /2; 0% = fEE 
For -u 34 we have: 

M 3 4 = £ + &3 4 X2 + (-/i,X2 + h- X i)E + bl 4 - 

For coefficients: 

684 = A & = /2; c& = fEE 
For -u 44 we have: 

« 44 = a?! - E z b° u + E+(- X 3Xi + X4X 2 )(a 44 + E z bf 4 ) + 
E + (-h z X3Xi + h'XiXi ~ h + X3X2 - Kx^X^ia^ ~ E z bf 4 ) + 
2XiX2X3X*(E + ) 2 (4 4 - E z b\ 4 ) 

For coefficients: 

a 2 4 \ = -f2EE; bf 4 = flEE; a 44 = flEE; bf 4 = fEE; 

a\ 4 = DEDEf 1; b\ 4 = DEDEf 

We see latter that many series in our list are equivalent to each other. After such selection 
between similar ones we have only some different series. 
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8 Analytical representation for series 



Collecting the terms in series expansion for a and b coefficients to 12 order we obtain for 
example the following representation for: 

1_ b 2 + h 2 b 4 + h 2 b 2 + h 4 b 6 + h 2 b 4 + h 4 b 2 + h 6 b 8 + h 2 b 6 + h 4 b 4 + h% 2 + h 8 
^ ~ 3! + ~~ 5! + 7! + 9! + lTl + 

6 10 + tftf + h 4 b 6 + h% 4 + h S b 2 + h 10 ^ 



13! 

b 12 + /^IQ + fr4 & 8 + h%Q + fr8ft4 + frl0fr2 + fe 12 

15! 

ftl4 + ^12 + ^10 + ^8 + ^6 + ^10^4 + ^12^2 + ^14 



+ 



17! 

Let us show how to sum following infinite series from b for example. We have 

b 2 b 4 b 6 b 8 b 10 , b 2 b 4 b 6 b 8 b 10 b 12 

1+ 2! + * + ^ + M + + ^ (1+ 3! + 5T + 7l + ^ + Tll + T3! + } 

It is seen that the first series equal to 

b 2 b 4 b 6 b 8 b w 
cosh(b) = l + - + - + - + - + - + 



and the second series equal to 

b 2 b 4 b 6 b 8 b W b 12 

sinh(b) b = 1 + — + - T + -T + - H H r + 

y " 3! 5! 7! 9! 11! 13! 

For sum of two series we have the following expression cosh(6) + E z sinh(6)/6 
The main series for us is the following expansion: 

1 + h 2 b A + fr2ft2 + h 4 b 6 + fr2ft4 + fe 4ft2 + fr6 g + fe 2ft6 + fr4fr4 + fe6& 2 + fr8 

* ~ 5! + 7! + 9! + ill + 13! + 

&10 + ^8 + ^6 + ^4 + ^2 + ^10 

+ 



15! 

b X2 + ^ftlO + W + h% S + fr8 & 4 + h X % 2 + fe 12 

17! 

&14 + ^12 + ^10 + ^8 + ^6 + ^10^4 + ^12^2 + ^14 



+ 



19! 

Let us show how the summation of this series can be performed. We can make the 
summation of subpart of hole series: 
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ft 5 b 7 V> b u b 13 , ,„ , If- 

v + v. + u + w + W + - = amh{b) - b -v 

h2 brV b« b" l?_bS 
v 7! 9! 11! 13! ' v v ; 3! 5! ; 

h %. + 11! + 13! + •- ) = * (8mh(6) " h - 3! " 5! " 7! } 
Having this series representation we can rewrite expression for f 

1 b 2 + h 2 b 4 + h 2 b 2 + h 4 b 6 + h 2 b 4 + h 4 b 2 + /i 6 b 8 + h 2 b 6 + h 4 b 4 + h 6 b 2 + h 8 



5! 7! 9! 11! 13! 

&10 + /^S + ^6 + h 6 b 4 + W + ^10 



+ 



15! 

&12 + fr 2& 10 + fr4ft8 + h% 6 + fe 8 & 4 + frl0 & 2 + fe 12 

17! 



+ .... = 



-(smh(6) - b - -) + -/> 2 (smh(6) - 6 - - - -) + ^ 4 (smh(6) - b - - - - - -) + .. = 

/? 2 /? 4 ft 3 /? 2 /? 4 

sinh(6)/6 5 (l + - + - + ...) + {-b - -)/b\l + W + ¥ + -..) + 

, b\h 2 , h 2 h 4 x , /i 2 h 4 

( -5! + * + ¥ + } + ( ~7! V* + ¥ + ¥ + ~> + - 

It is seen that the series of the type 1 + p- + |r + describe geometric series and 

gives the following result: 

h 2 h 4 1 b 2 

1 + TT + TT + 



6 2 ' 6 4 i_g 6 2 -/* 2 

If we insert this result we obtain the representation for: 



I b 2 + h 2 b 4 + fe 2fr2 + fr4 ft6 + fr2fr4 + fe 4 & 2 + fr 6 g + fe 2fr6 + fr4fr4 + h% 2 + fr8 

5! + 7! + 9! + 11! + 13! + 

ftlO + /^8 + ^6 + ^4 + W + ^10 

15! + - 



^ , , ft 3 ,,,, & 2 , & 5 > 2 & 2 , & 7 > 4 & 2 



b 2 -h 2 v 3! 7/ 6 2 -/i 2 v 5! 7 6 7 6 2 -/i 2 v 7\ J b 9 b 2 - h 2 
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sinh(b) 1 1 u2 



b 3 b 2 



1 h 2 h 4 , 1 



( Ql 71 "") 



6 2 -/i 2 b 2 -h 2 x 3! 5! 7! "' v 6 2 -/i 2 

sinh(6) 1 -. -. 



+ (-sinh(/i)//i' 3 + 
1 



&2_fc2 h 2_ h 2 ' y V // h 2J b 2 ~h 2 

sinh(fe) sinh(/i) 



6 3 _|_ 



Let us consider two series: one is / and second is f\. If we multiply / by coefficient 
a 5 and make following substitusion b— > ba, h— > ha we obtain the following series 

5 1 7 b 2 + h 2 q b 4 + h 2 b 2 + h 4 u b 6 + h 2 b 4 + h 4 b 2 + h 6 ^b 8 + h 2 b 6 + h 4 b 4 + h% 2 + h 8 
a V +a — +a 9! W 111 + ° 13! + 

ftlO + h 2 b S + ^6 + h% 4 + W + ^10 

° 15! + 

b 12 + ^10 + ^8 + h% 6 + ^4 + ^10^2 + ^12 

° 17! + 

It is obvious that if we take derivative we can reduce factorial in our series for example: 

_ d(a 5 f) ,dh , = dh = dh 

Jl o 1 J2 ; J3 1 J4 r% 

aa aa oa oa 

here we must insert b— > ab, h— > ah and a put to 1 after differentiation. 
Taking derivatives we obtain the analytical expression for f\ : 

^ co^) _ cosh^ ^ sinh^ ^ cOsh(&)-COsh(/l) _ 

h- b 2 -h 2 '> f 2 ~ b 2 -h 2 ] b 2 -h 2 ' 

6sinh(6) — hsinh(h) 
? 4 = W^h 2 
Series for have the following forms: 

1_ b 2 + h 2 b 4 + h 2 b 2 + h 4 b 6 + h 2 b 4 + h 4 b 2 + fr 6 6 8 + h 2 b 6 + W + h% 2 + fr 8 
^ ~ 4! + 6! + 8! + 10! + 12! + 

& 10 + h 2 b 8 + ^6 + h 6 b 4 + ^2 + ^10 

+ 



14! 

6 12 + tftfO + ^8 + h% e + ^4 + ^10^2 + ^12 



16! 

&14 + fr2frl2 + fr4ftl0 + h% 8 + h 8 b 6 + ^10^4 + ^12^2 + fe 14 

17! 



f 

i 

+ .... 



£ fe 2 + fr 2 fr4 + fe 2 & 2 + fr4 b 6 + fr2ft4 + fe 4ft2 + & 8 + ^6 + fr4 & 4 + h% 2 + fr8 

^ ~ 3! + 5! + 7! + 9! + 11! + 
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fr 10 + h 2 b 8 + h 4 b 6 + h% 4 + h 8 b 2 + h 10 

13! + 

b 12 + fr2 6 10 + fr4ft8 + fr6ft6 + h S b 4 + ^10^2 + fe 12 

15! + 

ftl4 + h 2 b U + ^4 fe 10 + ^8 + ^6 + ^10 fe 4 + ^12^2 + ^14 
1* + 

I b 2 + h 2 b 4 + h 2 b 2 + fr 4 b 6 + W + W + h 6 b 8 + h 2 b 6 + h 4 b 4 + h% 2 + h 8 
h ~ 2! + 4! + 6! + 8! + 10! + 

ftlO + ^8 + ^6 + ^4 + ^2 + ^10 

12! + 

6 12 + /^lO + ^8 + W + ^4 + ^10^2 + ^12 

14! + 

&14 + fe 2 & 12 + ^frlO + fr6ft8 + tftf + frl0fr4 + frl2ft2 + fe 14 

16! + 
U = l + — + 4! + 6! + 8! + 

6 10 + h 2 b 8 + h 4 b 6 + h 6 b 4 + h 8 b 2 + ^10 

10! + 

b 12 + tftfO + fr4ft8 + h%Q + fr8ft4 + ^10^2 + fr!2 

12! + 
b 14 + h 2 b 12 + h 4 b 10 + h% 8 + h 8 W + h 10 b 4 + h 12 b 2 + h 14 

14! + 

If we take the series for / and multiply it by 6 2 and take following derivative we 
obtain the series for f'EE. 

rrnrn 1 2b 2 + h 2 3b 4 + 2h 2 b 2 + h 4 4b 6 + 3h 2 b 4 + 2h 4 b 2 + h 6 

fEE =V + — — + 9! + U! + 

5b 8 + 4h 2 b 6 + 3h 4 b 4 + 2h% 2 + h 8 
13! + 
66 10 + 5h 2 b 8 + 4h 4 b 6 + 3h% 4 + 2h 8 b 2 + h 10 
15! 

For series fhh we must make multiplication of / on h 2 and make derivative on h 2 : 
fhh - 

1 2h 2 + b 2 3h 4 + 2b 2 h 2 + b 4 4h 6 + 3b 2 h 4 + 2b 4 h 2 + b 6 
/M= 5! + ^^ + 9! + 11! + 
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5h 8 + 4b 2 h 6 + 3b 4 h 4 + 2b 6 h 2 + b 8 
13! + 
6h 10 + 5b 2 h 8 + 4b 4 h 6 + 3b 6 h 4 + 2b 8 h 2 + b 10 
15! 

x ^ ^ 1 lb 2 + h 2 3b 4 + 2h 2 b 2 + h 4 46 6 + 3h 2 b 4 + 2h 4 b 2 + h 6 

jo hi hi = h h H h 

J 3! 5! 7! 9! 

5b 8 + 4/i 2 6 6 + 3/i 4 6 4 + 2h% 2 + /i 8 

+ 



11! 

66 10 + 5h 2 b 8 + 4h 4 b 6 + 3h% 4 + 2h 8 b 2 + h 10 



13! 

To evaluate fiEE, foEE^ fihh we multiply fEE by a 5 and make substitution 6— 
ba, h— > ha. After calculation we fix a = 1 

It is seen that expression for f\EE, faEE, f-Jih are the following: f\EE = ( %f g ) a =i 

r 1 2b 2 + h 2 3b 4 + 2h 2 b 2 + h 4 46 6 + 3h 2 b 4 + 2h 4 b 2 + h 6 

flEE = Ai + ^ + si + m + 

5b 8 + Ah 2 b 6 + 3h 4 b 4 + 2h% 2 + h 8 

+ 



12! 

6b 10 + 5h 2 b 8 + 4h 4 b 6 + 3h% 4 + 2h 8 b 2 + h 10 



+ 



14! 

Repeating all operation we can obtain for fiEE = ( 9 Q^ E ) a =i 

„ ^ ^ 1 2b 2 + h 2 3b 4 + 2h 2 b 2 + h 4 46 6 + 3h 2 b 4 + 2h 4 b 2 + h 6 

foEE = 1 1 1 h 

J2 3! 5! 7! 9! 

56 8 + Ah 2 b 6 + 3h A b 4 + 2h 6 b 2 + ^8 

+ 



11! 

66 10 + 5h 2 b 8 + 4h 4 b 6 + 3h% 4 + 2h 8 b 2 + h w 



and forflhh 



13! 

'dfhh^ 
v da 



la=l 



1 2h 2 + b 2 3h 4 + 2b 2 h 2 + b 4 4/i 6 + 3b 2 h 4 + 2b 4 h 2 + b 6 

flhh =A\ + — + 8! + 10! + 

5h 8 + 46 2 /i 6 + 3b 4 h 4 + 2& 6 fr 2 + & 8 | 



12! 

6/i 10 + 56 2 /i 8 + 4b 4 h 6 + 3b 6 h 4 + 2b 8 h 2 + b 10 



14! 



+ 



Series for DEDEf equal to derivative of fEE on 6 2 . It is seen if we compare series 
fEE and DEDEf: DEDEf = 

2 6b 2 + 2h 2 12b 4 + 6h 2 b 2 + 2h 4 20b 6 + \2h 2 b 4 + 6h 4 b 2 + 2/i 6 
^ /= 7J + ^^ + lii + 13! + 



16 



30& 8 + 20/i 2 6 6 + 12h 4 b 4 + Qh 6 b 2 + 2h 8 ^ 



15! 

We can obtain the following representation for DEDEfl : DEDEfl = 

2 6b 2 + 2h 2 126 4 + 6h 2 b 2 + 2h 4 206 6 + 12h 2 b 4 + 6h 4 b 2 + 2h 6 
DEDEfl = V + ^^ + 10! + 12! + 

30b 8 + 20h 2 b 6 + \2h 4 b 4 + 6h% 2 + 2h 8 



14! 

Comparing series for DEDhf2 and series for f 2 we can obtain DEDhf2 = B(h2 dl 



d(b 2 )d(h 2 ) 



, 1 2b 2 + 2h 2 3b 4 + Ah 2 b 2 + 3h 4 Ab 6 + Qh 2 b 4 + Qh 4 b 2 + Ah 6 

DEDhf2 = — H ■ 1 ■ 1 ■ h 

J 7! 9! 11! 13! 

5b s + 8h 2 b 6 + 9/i 4 6 4 + 8h% 2 + 5h 8 

+ 



15! 

66 10 + 10h 2 b 8 + 12h 4 b 6 + I2h% 4 + 10h 8 b 2 + Qh 10 



+ 



17! 

and for DEDhf3 we can obtain the following representation: DEDhf3 — — 



d(fc 2 )d(fe 2 ) 



„ „ , 1 26 2 + 2/i 2 36 4 + 4/i 2 6 2 + 3h 4 Ab 6 + 6h 2 b 4 + Qh 4 b 2 + Ah 6 
DEDhf3 =& + ^^ + lb! + 12! + 

56 8 + 8h 2 b 6 + 9h 4 b 4 + 8h% 2 + hh 8 | 



14! 

66 10 + I0h 2 b 8 + 12/i 4 6 6 + I2h% 4 + I0h 8 b 2 + 6/i 10 



+ 



16! 

Summing our calculation let's write all analytical formula for the function in the matrix 
elements: 

2 sinh(fr)6 3 + h (b 2 - h 2 ) cosh(6)6 + {h 3 - 3b 2 h) sinh(6) _ 
~ 2tfh{b 2 -h 2 ) 2 ' 

-26cosh(6) + 26cosh(/i) + (6 2 - /i 2 )sinh(6) 
~ 2b (b 2 - h 2 ) 2 ' 

, n ^ ^ b(b 2 - h 2 ) cosh(6) - [b 2 + h 2 ) sinh(6) + 2bh sinh(/i) 

j2EE = T5 

2b(b 2 -h 2 ) 2 

hsmh(b) — 6sinh(/i) cosh(fr) — cosh(/i) 

7 = ¥h^bh? ' 7 fe 2 "^! 2 ' 

6sinh(6)-fesmh(fe) 
7 6 2 -/i 2 ' 7 6 2 -/i 2 
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9 Analytical representation for supercoherent state 



Having nonlinear representation of underlying symmetry group of Hubbard model we can 
construct the supercoherent state as the action of the supergroup on some fix vector, for 
example | >- - the hight weight vector of the fundamental representation. 



G y= exp 



>- 



/ E z E+ \ 

xi K h + o 

X2 h~ -h z 
\ E~ xs Xa ~E z ) 

here | y is the eigenfunction of the Hubbard repulsion, describing the state with no 
electron. 

We write expression for | G >- in the following form: 



G y-- 



l 9o + 9ijXiXj + 94X1X2X3X4 \ 

gtxi + 92X2 + gt 2 4XiX2X4 

9lXl + 92X2 + 5 f l~23XlX2X3 

\ 9o+ 9ijXiXj + #4X1X2X3X4 / 



here coefficients g are equal to 



g% = cosh(fr) + E z sinh(6)/6; 



d 



sinh(6) sinh(ft) 



<?3°i = ~E + {a\\ + E z b\\) - E + h z (aH + E z b\\) = ~£ + (^y(& 2 ( 6 fi2 _ fe2 h )) + 



d 



sinh(6) sinh(?i) 



d 



sinh(6) sinh(/i) 



-)); 



g° 42 = E + (al\ + E z b\\) - E + h z (aH + E z bH) = E~ 



d 



'd(b 2 ) 



sinh(6) sinh(ft) 

b 2 _ h 2 )n 



o sinh(b) sinh(?i) 



d 



sinh(6) sinh(?i) 



-)); 



d 



sinh(6) sinh(ft) 



g° A1 = E+h-(a? 1 + E z bH) = ^"(^(^(^ ^ 



o sinh(6) sinh(ft) 



rs sinh(6) sinh(ft) 

g° 32 = -E+h + (ali+E z bH) = E+h + {J-Ab^ 



dsinh(6) sinh(ft) 



d{b 2 ) 



b 2 ~h 2 



d(b 2 



b 2 -h 2 



))); 
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o2 sinh(b) sinh(fe) rj2 sinh(6) sinh(/i) 

a sinh(6) sinh(/i) o sinh(6) sinh(/i) 

ffl =/i (a 31 + £*6 31 ) = /i ( &2 _ fe2 + &2 _ fe2 ), 

/ n ^ ,n n , / i ^ ,1 n 6sinh(6) — /isinh(/i) 
g 2 = (a° 31 + f?^) - h z {a\i + #*&3i) = ~jrZ^ ~+ 



b 2_ h 2 ~* * b 2_ h 2 ' 

o sinh(6) sinh(fe) n sinh(6) sinh(ft) 

gl = E~ sinh(6)/6; 



.2 2.1 ,-+,-,,21, , ,,.22 L - +L -,22, _ cosh ( & ) _ cosh (/ l ) 



<& = -« + E + E~bi[) - h z (aH + E + E-bji) = -h z 



b 2_ h 2 > b 2_ h 2 



o cosh(6) cosh(/i) o sinh(b) sinh(/i) 

o cosh(fe) cosh(ft) sinh(6) sinh(ft) 
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sinh(b) sinh(ft) o sinh(6) sinh(ft) 

& = h-(a% + E+E-%) = M ~V^ ) + ^ + ^"^) I 



sinh(b) sinh(ft) o sinh(6) sinh(ft) 

gl 2 = -h+(all + ^S-6g) = - ^ + ^ + ^y(^ 2 ( ^:^ )); 

o sinh(fe) sinh(ft) 

, 4 2 = 2£+(c4 + E+E-biJ = 2^ + (^y(5 2 ( ^:^ )) + 

02 sinh(b) sinh(/i) 

10 Conclusion 

We have calculated the exact representation of the supergroup as well as the supercoherent 
state in the Hubbard model. These constructions naturally appear in the strongly correlated 
electronic systems in the case of introducing atomic base in limit of large on-site Hubbard 
repulsion. The dynamical supergroup which operates in a local superbundle determined by 
any on-site eigenfunction gives us the wave function in the form of a superspinor. This super- 
spinor describes a local supercoordinates frame in the curved supermanifold. The operator 
spinor part acting in tangent and cotangent bandies of this supermanifold in supergroup 
can be reformulated in the terms of the atomic Hubbard operators. Next step lies in the 
calculation of effective functional of Hubbard model. 
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